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Abstract 

Non-generic contributions to the quantal level- density from parallel segments 
in billiards are investigated. These contributions are due to the existence of 
marginally stable families of periodic orbits, which are structurally unsta- 
ble, in the sense that small perturbations, such as a slight tilt of one of the 
segments, destroy them completely. We investigate the effects of such per- 
turbation on the corresponding quantum spectra, and demonstrate them for 
the stadium billiard. 
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1 Introduction 



Hamiltonian systems which display strong forms of chaotic dynamics (K- 
systems) may tolerate families of periodic orbits which are marginally stable. 
They occupy a domain of measure zero in phase space and therefore do not 
affect the ergodicity of the system as a whole. These measure zero sets of 
periodic orbits are non-generic and also may be "structurally unstable" in 
the sense that an infinitesimal perturbation may destroy their periodicity 
completely. This sensitivity is in stark contrast with the structural stability 
of the chaotic component, which is robust against small perturbations. A 
notable example of such sets of trajectories are the "bouncing ball" (BB) 
orbits in the stadium billiard [1], which bounce perpendicularly between the 
parallel straight sections of the boundary. It is a continuous one parameter 
family of neutral periodic orbits, which disappear when the straight segments 
are slightly tilted. 

Even though the non-generic components play insignificant role in the 
classical description, they might have rather noticeable effects on the spec- 
tral properties of the corresponding quantum analogues. Berry [2] studied 
their effect on the spectrum of the Sinai billiard. More recently they were 
brought again into the focus because of their prominent appearance in the 
spectrum of the stadium billiard [3, 4]: the Fourier transform of the spectral 
density (obtained numerically [4] and experimentally [3]) shows a sequence 
of structures which correspond to the multiple traversals of the bouncing ball 
orbits. This phenomenon was analyzed and explained in subsequent papers 
by Sieber et al [5] and by Alonso and Gaspard [6]. In view of the structural 
instability of these orbits, one could naively expect that the corresponding 
quantum features would also disappear as soon as the boundary is slightly 
deformed. However, quantum dynamics is a wave phenomenon, and as such 
it cannot resolve deformations which are smaller than the shortest wave- 
length in the relevant domain. Thus, it is more reasonable to assume that 
the quantum features which are due to non-generic and structurally unstable 
classical orbits will disappear gradually when the system is perturbed, and 
will leave their strongest mark on the higher parts of the quantum spectrum. 

The theory in the next section, as well as the computational results of 
section 3, will explicitly demonstrate that the above expectation of gradual 
rather than abrupt transition actually holds. Also, a quantitative measure 
for the amount of deformation that is needed in order to observe the change 
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in the quantum spectrum will be derived. 



2 Theory 

We shall start the discussion by deriving a general expression for the semi- 
classical contribution to the density of states function, produced by classical 
orbits that bounce from straight components of a planar billiard, exclusively. 
We shall apply this expression for billiards with parallel straight segments, 
such as e.g. the Bunimovich stadium [5], and the Sinai billiard [2]. The 
emphasis will be put on identifying the contributions that are generated by 
the classical periodic orbits. Then, the classical dynamics of orbits that 
bounce between two slightly tilted segments will be analyzed, and the results 
will serve to derive the resulting non-generic contribution to the density of 
states. 



2.1 Level Density 

The standard Gutzwiller's trace formula [7], which is frequently used to cal- 
culate the level density semiclassically for chaotic systems, does not account 
for the special effects of non-generic periodic orbits [5]. Therefore, we start 
our discussion from the expression for the more fundamental semiclassical 
propagator. 

The contribution to the semiclassical propagator due to classical paths 
which reflect only from straight segments of the boundaries is given by: 
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where m is the mass of the particle, L is the length of the classical path from 
q 1 to q ", and v is a phase index, which is twice the number of collisions 
with the boundaries (for Dirichlet boundary conditions). It is assumed that 
between successive bounces, the motion is free. 

The propagator can be Fourier-transformed to give the outgoing Green's 
function: 

G"{q',q",E) = lim [°° K"(q ',q ",t)exp I Ue + ie) t\ dt (2) 
e^o+ in Jo [n J 
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" m 2 E ^ lj (^(9',9"))exp{-z-^ (3) 

where Hq 1 ^ is the Hankel function of order zero and first kind [8], and k = 
\l2raE /h is the wavenumber. Note, that for this special case, the integral 
in (2) could be calculated exactly, without invoking the stationary phase 
approximation. 

The level density function is defined by: 

d(E) = --5s [ G(q,q,E)dq (4) 

7T J 

and we define for convenience the function d(k) } such that d(k)dk = d(E)dE. 
Inserting eq. (3) into (4) we get: 



d s (k) = —d s (E) = ^f J2 Mkmd))eM-^W (5) 
m Z7T J , Z 
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where the summation extends over all classical paths from qtoq that reflect 
from straight segments only. Since d(k) is a linear functional of the propa- 
gator K, one can identify d s (k) as the semiclassical contribution to the level 
density of those orbits that reflect from straight segments only. 



2.2 Application to the Two-Parallel Lines Case 

As a first application of (5), we shall study the case of two parallel lines of 
length b and at a distance a, such as can be found in the stadium [1] (see 
fig. 1). The results of this calculation will serve as a reference point for the 
more general case to follow. In order to perform the calculation of d s (k) } we 
should find for every point q inside the domain all the classical orbits that 
begin and end in q. It is not difficult to see that these orbits are the ones that 
impinge vertically on the parallel lines. Actually, there is an infinite number 
of such orbits for every point q, and we can divide them into two different 
classes: 

1. Periodic or "even" orbits — these are the orbits whose final momenta 
are parallel to their initial momenta. They bounce an even number of 
times, and therefore have an overall phase of +1. 
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2. "Odd" orbits — these are the orbits whose final momenta are anti- 
parallel to their initial momenta. They bounce an odd number of times, 
and their phase is -1. 

The corresponding lengths of these orbits are: 

l/2 

L e (n,x,y) = [(x 2 - x x f + (y 2 - yi + 2cm) 2 ] , n £ Z (6) 

l/2 

L (n,x,y) = [(x 2 - x x f + (y 2 + yi + 2cm) 2 ] , n G Z. (7) 

where we have explicitly used the Cartesian coordinates: q = (x, y). Inserting 
the above lengths into equation (5) we get after some manipulations: 

abk + °° b 
d s {k) paraUe i = — V J {2akn) - — = d e (k) + d (k) (8) 

Z7T — Z7T 

n= — 00 

This result has been previously obtained by Berry [2], Sieber et al [5] and 
Alonso and Gaspard [6]. It differs from the generic contributions to d(k) 
that come from unstable, isolated periodic orbits in two respects: firstly, it 
is composed of contributions from continuous families of neutral trajectories. 
Secondly, it includes contributions from orbits that do not close in phase 
space, but rather close in configuration (coordinate) space. 

Analyzing eq. (8), we notice that the periodic (even) orbits give oscilla- 
tory terms in k. Actually, if d e (k) is Fourier-transformed [5], strong peaks are 
found for lengths ±2ma (m is a positive integer), corresponding to the clas- 
sical periodic orbits . These peaks have also been observed in the microwave- 
cavity experiment [3] and in numerical simulations [4]. The n = term in 
(8), which is the contribution of zero-length orbits, gives a smooth (non- 
oscillatory) contribution abk/2TT. This is the leading (area) term in Weyl's 
asymptotic expression for the smooth level density in planar billiards [2, 9, 5] 
coming from the rectangle a X b. The odd orbits yield a constant term, which 
is the contribution of the parallel lines to the first correction (perimeter) term 
in Weyl's formula. 

There is a difference in powers of k (and thus in powers of fi) between the 
contributions of the even and the odd orbits in (8). It is due to interference 
between the contributing amplitudes, which is constructive for the even or- 
bits, and is destructive for the odd orbits. In the latter case there is an exact 
cancelation of the terms of order \fk. 
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2.3 Classical Dynamics of the Deformed Billiard 



The classical trajectories which bounce between the parallel segments are 
structurally unstable. An infinitesimally small tilt of the lines away from 
perfect parallelism, will cause the continuous family of periodic orbits to 
disappear. This, in turn, means, that had we used Gutzwiller's formula 
[7] (or some other variant that relies upon periodic-orbit theory) to calculate 
d(k), a discontinuous behavior (in some displacement parameter) would have 
shown up. In the following we shall show that the transition due to such a 
tilt is actually a continuous one. 

We shall follow the behavior of the non-generic expression (8) when the 
straight segments are tilted at a small angle <p (fig. 2). This justifies a first- 
order analysis, which makes the calculations tractable. 

Once we have tilted the lines in the way described, the classical orbits 
which go from q to q will be deformed too: while in the parallel case the 
x-coordinate of the contributing orbits remains the same for all points of 
the orbit, it will be no longer the case in the deformed billiard, and the 
trajectories will stray from their initial x-coordinate. 

Consider all the trajectories which start from an arbitrary point between 
the lines. It is easy to see, that if the trajectory initially aims at the narrowing 
part of the billiard, the scattering angle between the straight segment of the 
boundary and the trajectory will gain an angle <p at every reflection from the 
lines. Therefore, the total scattering angle will eventually exceed tt/2 and 
the trajectory will turn around towards the widening part of the billiard. 
This suggests, that there might be some initial launching angles, for which 
the trajectory will come back to its initial point, and a closed trajectory will 
be traced. If, however, the trajectory will initially go to the widening part, 
it will escape without ever coming back to its initial position. 

We shall use elementary geometry to identify the angles for which the 
trajectories close. We shall first consider trajectories with initial momentum 
in the "upward" (negative y) direction. Since the scattering angle increases 
by (p for every collision with the straight segments, it is not difficult to see, 
that a trajectory which starts with an angle of rrup (where m is a non-negative 
integer) relative to the negative y axis, will scatter from the boundary at an 
angle tt/2 after m + 1 collisions. Hence, the trajectory will retrace itself, 
coming back to its initial position with final momentum that is opposite to 
the initial one. Fig. 3 shows such a trajectory, for m = 1. Thus, the launching 
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angles: 



a n 0+ = nip , n = 0, 1,2, . . . 



(9) 



yield closed orbits. The subscript o+ stands for the fact, that these orbits 
substitute the odd orbits that start upwards in the parallel lines case. We 
note that the selected values of the launching angles are independent of the 
position <f, but for every point there is an upper limit for the launching angle 
(due to the finite extent of the straight segments). Also, no approximation 
had been used in obtaining this set. 

In addition to the above family of closed orbits, we can find another 
(non-retracing) family that starts upwards (see fig. 3). Simple geometric and 
continuity arguments show, that between every two orbits with launching 
angles of rrup and (m + l)ip, there is a unique closed orbit that closes from 
below. That is, its final momentum is approximately parallel (pi ■ pj > 0) 
to the initial momentum. First-order (in the tilt angle) analysis gives the 
following formula for the launching angles: 



where the subscript e stands for the fact, that these orbits substitute the 
even orbits in the parallel-lines case. Here the launching angle depends on 
the y coordinate. 

Until now we have considered only trajectories with initial momentum in 
the upward direction. If we include now trajectories with downward initial 
momentum, we get an additional family of odd-type trajectories: 



where the launching angle is now being measured relative to the positive y 
direction. The absence of the n = term in d a _ reflects the asymmetry of 
the coordinate system relative to the inclined lines, and does not have any 
special physical meaning. As regard the other family of orbits (eq. (10)), the 
uniqueness of the trajectories and forward-backward time symmetry implies 
that no new trajectories of this kind can be found. The same geometrical 
paths should be used, with the direction of propagation reversed. 




, n = 1,2, . . . 



(10) 
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a n _ = rap , n = 1,2, . . . . 
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First-order calculations in the tilt angle <p give the following expressions 
for the lengths of the closed orbits: 

L n 0+ = 2{y + n[a + <p(b-x)]} , n = 0,l,2,... (12) 
L n _ = 2{n[a + v(b-x)]-y} , n = l,2,3,... (13) 
L n e = 2n[a + <p(b-x)] , re = 0,1,2,... (14) 

These expression are valid as long as the maximum excursion of the closed 
orbit in x direction is small compared to b. 

2.4 Calculation of the Non-Generic Contributions to 
the Level Density in the Deformed Billiard 

We shall first study the contribution of the even family of orbits, and later 
will give some arguments, why the odd family can be ignored. 

The calculation of d e (k) requires the application of eq. (5) to the family 
of orbits given in (14). However, special care should be given to the limits of 
the domain of integration in (5). 

The relevant domain can be divided into 4 different parts, shown in fig. 2. 
Region 1 is the original a X b rectangle, and stands for the main contribu- 
tion to the integral. Region 2 is formed due to the additional area in the 
deformed billiard relative to the original (parallel) one. Region 3 stands for 
the points in region 1, from which the classical trajectories eventually hit 
the chaotic scatterer. We should eliminate these points from the integration, 
since they no longer belong to the continuous family of closed orbits. Region 
4 is composed of all the points that reside outside region 1, which are starting 
points of closed orbits which scatter only from the straight segments. The 
boundaries of regions 1 and 2 are independent of the number of collisions re 
and the family of trajectories. The boundaries of regions 3 and 4 depend on 
both re and the family of trajectories. They can be calculated (to first order 
in (,c), and the result for the even family is: 

{ b + <p(£-y-*£) re even 
€L(n,y) = \ (15) 
( b + ip(£- ^±11) n odd 
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X 



(4) 

even 



2n-l 



[n,y) = { (16) 
<p(™ - (n - l)y - y > ^ 

'ELi^IbL are the boundary curves for regions 3, 4, respectively. 
We consider first the contributions from regions 1 and 2. The contribution 
of the even orbits is calculated by inserting (14) into (5). Recalling the 
two families of even orbits, and the +1 phase generated by even number of 
bounces, it will be given by: 



d e (k) = 2Y j ^- t dx r + ' Ab X) dyj {2nk[a + p(b-x)]} (17) 

n=1 ZTT Jx=0 Jy=0 



where N e stands for the number of valid terms to be included in this first- 
order analysis. Note that the domain of integration includes regions 1 and 2 
only. The y integration is trivial, and introduced a factor of a-\-p{b— x) (which 
is the local width) to the integrand. Changing variables to t = a + pib — x), 
the integral can be performed analytically: 

b 1 (a + bp) ^ [2nk(a + bp)] - aJ^nka) 

d e{k) = —l^ 7 • 18 

ZTT ^ n o<p 

Each term in the sum is in the form of a "finite derivative" of the function 
tJ\(t) with respect to b<p. Taking the limit — > and using the identity 
tJ (t) = tJ[(t) + J\(t) } give us exactly the first term in (8), as expected. 

The contributions from regions 3 and 4 are linear in (p. Since the main 
effect is dephasing between two leading-order terms (18) we can neglect the 
contributions of regions 3 and 4. This was validated by analytic and numer- 
ical calculations. 

A better insight to the properties of d e (k) is gained by considering the 
length spectrum: 



'"max 



D e (x, k max ) = — — I cos(kx)d e (k)dk. (19) 

The domain of integration is taken to be finite (after [5]), in order to make 
it possible to relate the results to experiments. Inserting (18) into the above 
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definition (19), and using the asymptotic expression for Ji(x) } x > [8] we 
get: 



d e {x , k max ^ — ^ 



\ 4:Trk max n(fi^/2Trn 

S ' (k max A^~) C(A; ma:E A~'~) S(k max A-) + C(k 



1 a + hip 



S{k m axB~^) — C (k max B^~) S{k max B ) + C {k max B ) 



where we have defined: 



A ± 
B ± 



B+\ 



x ± 2n(a + b(p) 
x ± 2na 



B- 



d: 



S^) 



sign(xj 



M sin(t) 



o 

M cos(t) 



dt 



o 



dt. 



Prominent peaks in (20) are expected for vanishing denominators: x = 
±2n(a + b(p) , x = ±2na. Treating the x > domain, we have the fol- 
lowing picture. For very small (p } the two peaks are very large (of order 
l/(fi) and very close (of order <p) to each other. Since D e in (20) has the 
structure of a "derivative", we expect to observe two sharp extrema with 
opposite signs near x = 2na } much like the case of differentiating a Gaussian 
or a Lorentzian. When <p increases, the two peaks separate, their amplitude 
decrease, and we obtain destructive interference between the terms when: 



2nk max bcp tt. 
This defines the critical tilt angle: 

TC 



7T 



2Tik max b 



(21) 



(22) 



20) 
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for which we expect the structure of D e near x = 2na to be appreciably lower 
and broader than for the <p — > case. 

A similar analysis had been applied to the odd family of orbits. The main 
result was, that D (x } k max ) is sharply concentrated near x = for wide range 
of tilt angles, with no other prominent peaks. Since there are other contri- 
butions which are independent of k (e.g., Weyl's term from the perimeter), 
it can be difficult and not practical to resolve the specific contribution of the 
odd family from all the other contributions. 



3 Numerical Results and Discussion 

In this section we shall demonstrate numerically the behavior of D e for a set 
of parameters that had been used by Graf et al [3] to measure 1060 levels 
of the Bunimovich stadium, using superconducting microwave cavity. (The 
actual experiment was carried out in a parallel geometry, and we use the 
same parameters just for convenience.) The values of the parameters are: 
a = 20 cm, b = 36 cm, k max = 3.665 cm -1 . Substituting these values into 
eq. (22) for the critical tilt angle we get: 

« (23) 

n 

In fig. 4 we plot D e (x } k max ) (calculated from (20)) for three different 
values of (p } and for the term n = 1. For an extremely small angle, <p = 
10~ 6 <C fl r , we recover the results of Sieber et al. [5] for the parallel case, 
and the function has two prominent extrema with opposite signs near x = 2a 
= 40 cm, as anticipated above. For <p = 0.012 = tp\ r some lowering and 
broadening is observable. This trend is fully developed for <p = 0.03, where 
the peaks had been widen considerably, resulting in an amplitude which is 
1/3 of the original one. Fig. 5 shows the heights of the positive peaks of 
D e (x } k max ) as a function of <p for n = 1 and n = 2 in (20) (the values 
are normalized according to the <p = case). The decline of the maximal 
amplitude is evident, and the critical angle (23) (marked by cross) gives a 
good estimate for this decline to be already significant. We can also observe, 
that the decline for n = 2 is twice faster than for the n = 1 case. This is 
consistent with our prediction of the critical angle, eq. (23). 

In summary, we obtained in this work a quantitative measure of the scale 
of deformations which are necessary to eliminate the effects of the bouncing 
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ball orbits. This is of interest in the semiclassical study of billiards [5], where 
the non-generic families of orbits are a nuisance which masks and interferes 
with the contributions of the unstable periodic orbits. The effect can be 
demonstrated in microwave cavity experiments of the type reported in [3]. 
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Figures 

Fig. 1: Geometry of the stadium billiard. 

Fig. 2: Integration regions for the tilted stadium. 

Fig. 3: Closed orbits in the tilted stadium. Right: odd orbit, left: even orbit. 

Fig. 4: Length spectrum D e (x } k max ) for n = 1 and tilt angle <p = 10~ 6 , 0.012 and 0.03 rad. 

Fig. 5: Maximum of D e (x } k max ) as a function of tilt angle <p for n = 1 and n = 2. 
Values are normalized according to <p = case and crosses mark <p™ r . 
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Fig. 5 




19 



